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Introduction 


The  purpose  of  this  note  is  to  outline  a  set  of  procedures  for  processing 
hourly  observations  of  tidal  elevation  to  obtain  information  on  subtidal 
variabi 1 ity . 


Recommended  Processing  Procedures 

Sea  level  data  obtained  from  the  National  Ocean  Survey  of  the  National 
Oceanic  and  Atmospheric  Administration  (NOAA/NOS)  tide  gages  is  normally 
provided  on  an  hourly  basis  with  tidal  elevations  given  in  feet.  The  first 
correction  to  be  applied  adjusts  the  tidal  elevations  which  are  referenced  to 
an  arbitrary  datum  initially,  to  datums  which  are  referenced  to  either  Mean 
Higher  High  Water  (MHHW),  Mean  High  Water  (MHW),  Mean  Sea  Level  (MSL),  Mean 
Low  Water  (MLW),  or  Mean  Lower  Low  Water  (MLLW).  NOAA/NOS  provides  the 
necessary  data  for  each  tide  gauge  location  to  make  these  corrections.  The 
correction  takes  the  following  form: 

Adjusted  Sea  Level  (feet)  =  Raw  Observation  (feet)  -  Correction  (feet) 

To  obtain  MSL  for  example, 

MSL  (feet)  =  Raw  Observations  (feet)  -  MSL  Correction  (feet) 

It  is  noted  that  because  relative  variations  in  sea  level  are  generally  of 
primary  interest  with  regard  to  subtidal  variability,  no  formal  requirement  to 
apply  a  datum  correction  in  this  case,  exists. 

The  second  adjustment  that  must  be  applied  removes  the  diurnal  and 
semidiurnal  tidal  components.  The  Godin  low-pass  digital  filter  is  used  to 
remove  these  components  (Godin,  1972).  The  response  of  this  filter  can  be 
expressed  as 
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where 


'2 

A  •  A  /25-242 
25  24 


A  =  sin(mTAta) 

^  '  s  1  n(iTAaa  ) 

0  =  frequency  in  cycles  per  day 
At  =  time  interval  between  observations  in  days 
Ap  represents  the  arithmetric  summation  of  n  observations  with  the  filter 
operating  on  the  data  three  times,  first  using  25-point  means  )  and  then 
using  24-point  means  twice.  According  to  Godin,  the  frequency  band  of 

this  filter  extends  from  0  to  0.8  cycles  per  day  (or  30  hours),  beyond  which 
the  response  is  effectively  zero.  The  response  of  this  filter  is  shown  in 
Figure  1.  The  half -ampl itude  cutoff  is  approximately  2.8  days  (or  -66  hours). 
The  amplitude  of  the  diurnal  and  semidiurnal  tidal  components  is  reduced  to 
less  than  2%  of  the  maximum  value  within  the  passband.  A  FORTRAN  program 
implementing  this  tide-suppressing  filter  is  contained  in  Appendix  1. 

A  third  adjustment  must  be  made  to  sea  level  data  to  remove  the  effects 
of  atmospheric  pressure.  According  to  Lisitzin  (1974),  in  the  stationary 
case,  the  sea  surface  is  depressed  by  1  cm  for  every  increase  in  atmospneric 
pressure  of  1.005  mbar  (and  vice  versa).  It  is  noted,  however,  that  such  a 
simple  static  adjustment  for  atmospheric  pressure  does  not  take  into  account 
the  effects  of  surface  wind. 

A  relationship  of  the  following  form  is  recommended  for  making  the  sea 
level  pressure  (SLR)  correction: 


SL(t)artj  =  (SL(t)-'Sr)  -  {SLP(t)-TOF) 


(1) 


*  1.000  vice  1. OOOcm/1.005  mbar  is  used  here 
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where 


SL(t)acij  =  Pressure  corrected  sea  level  (cm) 

SL(t)  =  filtered  sea  level  (cm) 

SL  =  mean  value  of  filtered  sea  level  (cm) 

SLP(t)  =  filtered  sea  level  pressure  (mbar) 

SLP  =  mean  value  of  filtered  sea  level  pressure  (mbar) 

The  filtered  sea  level  should  first  be  converted  to  centimeters.  Next, 
the  SLP  data  should  be  low-pass  filtered  in  a  manner  similar  to  that  used  for 
the  sea  level  data.  In  this  regard,  it  is  recommended  that  the  Godin  filter 
be  applied  to  the  SLP  data  as  well.  Equivalent  filtering  of  SLP  has  the 
advantage  of  producing  a  smooth  time  series  with  the  same  "time  constant"  or 
degree  of  smoothness  as  the  sea  level  data.  As  indicated,  both  SL(t)  and 
SLP(t)  should  be  mean  corrected  in  order  to  complete  the  correction  for 
atmospheric  pressure.  Finally,  whether  or  not  the  two  terms  on  the  R.H.S.  of 
(1)  are  actually  added  or  subtracted  will  depend  on  the  sign  convention  used 
to  define  sea  level  and  sea  level  pressure.  The  correct  choice,  of  course, 
will  produce  an  adjusted  series  with  decreases  (increases)  in  sea  level  that 
correspond  to  increases  (decreases)  in  sea  level  pressure. 

An  Example 

A  92-day  sequence  of  hourly  sea  level  data  from  the  NOAA/NOS  tide  gauge 
at  Santa  Monica,  California  was  acquired  to  illustrate  the  foregoing 
procedures.  A  datum  correction  furnished  with  the  data  from  NOAA/NOS  was 
applied  initially  to  produce  a  sequence  referenced  to  MSL.  Next,  the  Godin 
low-pass  filter  was  applied  to  the  sea  level  data  to  remove  the  diurnal  and 
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semidiurnal  tidal  components  (Figure  2).  The  upper  panel  in  Figure  2  shows 
the  mixed  character  of  this  tidal  record  which  clearly  contains  both  diurnal 
and  semidiurnal  components.  The  lower  panel  shows  the  sea  level  record  with 
the  tides  removed.  To  apply  the  correction  for  atmospheric  pressure,  hourly 
SLP  data  from  a  nearby  MDBC  environmental  data  buoy  (NDBC  Buoy  46011  off  Pt. 
Conception)  were  acquired.  These  data  were  likewise  smoothed  using  the  Godin 
low-pass  filter  (Figure  3).  The  predominant  diurnal  component  in  SLP  is 
effectively  removed  using  this  filter  (lower  trace).  In  the  final  figure 
(Figure  4),  the  mean-corrected  SLP  (upper  panel)  is  subtracted  from  the 
mean-corrected  sea  level  (middle  panel)  to  yield  sea  level  which  has  been 
corrected  for  atmospheric  pressure  (lower  panel  ). 

Discussion 

Many  different  factors  influence  sea  level.  The  time  scales  over  which 
these  factors  are  important  vary  from  seconds  to  millenia.  Only  relatively 
short-term  effects  on  sea  level  (<  a  few  days)  have  been  considered  in  the 
foregoing  discussion.  Although  use  of  the  Godin  filter  was  recommended  for 
removing  the  semidiurnal  and  diurnal  tides,  other  tide-removing  filters  exist 
(Walters  and  Heston,  1982).  In  fact,  Walters  and  Heston  found  that  a 
so-called  transform  filter  they  describe  generally  had  more  desirable 
properties  than  the  Godin  filter.  The  only  significant  drawback  they  found  in 
using  the  Godin  filter  is  that  in  addition  to  removing  tidal  variations,  it 
also  attenuates  variations  in  the  2  to  4  day  range.  For  applications  where 
information  in  this  particular  band  is  not  of  primary  interest,  however,  the 
Godin  filter  is  quite  satisfactory. 

The  simple  adjustment  for  atmospheric  pressure  described  here  assumes 
that  the  direct  response  of  the  ocean  to  a  change  in  SLP  is  hydrostatic. 


Hamon  (1966),  however,  has  shown  that  such  a  simple  one-to-minus  one 
correspondence  between  SLP  and  sea  level  may  not  always  be  observed.  Factors 
such  as  coastal  trapped  waves  may  contribute  to  the  deviations  from  the 
hydrostatic  approximation  (Hamon,  1976).  It  is  also  possible  to  examine  the 
phase  relationship  between  SLP  and  sea  level  using  the  techniques  of 
cross-spectral  analysis  (Hamon,  1966).  When  such  an  analysis  is  conducted, 
corrections  for  atmospheric  pressure  can  be  defined  for  a  range  of 
frequencies.  Using  this  approach,  Hamon  (1966)  found  that  changes  in  sea 
level  often  follow  changes  in  SLP  with  lags  ranging  from  0  to  2  days. 
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APPENDIX  1 
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GODIN  FILISF  A??LIES  A  24  2F  bCCCAF  I«IC"  AVD  A  2 
CNCE  TO  ■-■'C'JF.LI  DA2A.  T.-’S  FIPST  GOTP'n  FCINI  IS  3 
liiAN  Tils  pr-S''  IND’JT  DCINT. 


CIlENoION  51(3333) 

INPUT  ACTUAL  'lUXL-P 
SIATE^ENI 
N.TEAD  =  NU-iJES  C?  INPUT  DATA  tCINT-S 
DO  b6  J=1,NPSAD 
RSAC(5,11)  SL(J) 
bo  C'CNIINUS 

INSEFl  ACTUAL  TATA  PO’^.rAT  I'l  THE  PCLLC^ING  AT 
11  fCP.'l  AT  (D  AT  A  FCD.'IA:) 

FILTER  THE  DATA 

CALL  3CXCAr(SL,NFEAD,2i,1,NCCT1) 

CAIL  aC:<CAF  |SL,  NC'JTI  ,24,  1,  \CUT2) 

CALL  BCXCAR  ( S  L  ,  NC  UT  2 , 2  ,  1  ,  NC  U  T  3) 

T'Ie  FILISStD  l^7A 
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WRITE  OUT 
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CO  b7  i<  =  1,NCII'rj 

WRITE  (6, 12)  SL(<) 

b7  CONTINUE 
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END 
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33  CONTINUE 
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